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Weak magnetism increases antineutrino mean free paths in core collapse supernovae. The parity 
violating interference between axial and vector currents makes antineutrino-nucleon cross sections 
smaller then those for neutrinos. We calculate simple, exact correction factors to include recoil and 
weak magnetism in supernova simulations. Weak magnetism may significantly increase the neutrino 
energy flux. We calculate, in a diffusion approximation, an increase of order 15% in the total energy 
flux for temperatures near 10 MeV. This should raise the neutrino luminosity. Weak magnetism 
also changes the emitted spectrum of Vx (with x = ^ or t) and i'e- We estimate that Vx will be 
emitted about 7% hotter than because Vx have longer mean free paths. Likewise weak magnetism 
may increase the Ve temperature by of order 10%. This increase in temperature coupled with the 
increase in neutrino luminosity should increase the heating in the low density region outside of the 
neutrino sphere. This, in turn, could be important for the success of an explosion. It is important 
to check our results with a full simulation that includes Boltzmann neutrino transport and weak 
magnetism corrections. 



I. INTRODUCTION 



Core collapse supernovae are dominated by neutrinos, therfore many supernova properties may depend on the 
nature of neutrino- nucleon interactions. This provides an opportunity to study characteristic symmetries and features 
of the Standard Model weak interactions. Supernovae may provide macroscopic manifestations of charge conjugation 
||,| or parity violation 

In the laboratory at high energies, antincutrino-nuclcon cross sections are systematically smaller than neutrino- 
nucleon cross sections. This is related to charge conjugation and parity violation in the Standard Model. Nevertheless, 
most core collapse supernova simulations still use the same lowest order cross section dao/dQ for both neutrinos and 
antineutrinos, 

^ = ^[clil + COS0) + cli3 ^ cos9)]. (1) 

Here k is the (incoming) neutrino energy, 6 the scattering angle and the vector c„ and axial vector Ca coupling 
constants are listed in table I. 

In this paper, we discuss free space corrections to Eq. (^ from nucleon recoil, weak magnetism, strange quarks and 
single nucleon form factors. These corrections may be important for supernova simulations because they are present 
at all densities. In contrast, the density dependent nucleon correlations discussed by Q and [Q arc important only at 
very high densities, well inside the neutrino sphere. 

Electron antineutrinos were detected from Supernova 1987A via their capture on protons. Weak magnetism re- 
duces this cross section and impacts the deduced neutron star binding energy and antineutrino temperature. Weak 
magnetism also changes the rate of antineutrino capture in the neutrino driven wind above a protoneutron star. This 
may change the electron fraction Yj, (ratio of electrons or protons to baryons) and nucleosynthesis yields in the wind 
ii- 

However, to confirm this change in Ye, one should perform supernova simulations including weak magnetism because 
this may change the neutrino luminosities and emitted spectra which will also impact Y,,. In this paper, we present 
simple formulae so that these corrections can be incorporated in simulations. We also give estimates of the change in 
emitted electron antineutrino spectrum and the differences expected between the spectrum of mu or tau antineutrinos 
and neutrinos. 
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The reduction in opacities from recoil and weak magnetism should increase luminosities of both neutrinos and 
antineutrinos as more energy is transported to the neutrino sphere. This could increase heating behind the shock 
and more then compensate for the smaller antineutrino absorption cross section. Weak magnetism should also reduce 
the cooling from positron capture on neutrons further increasing the net heating behind the shock. These changes in 
heating could make a simulation more likely to explode. 

Finally, weak magnetism allows mu and tau antineutrinos to escape the star faster than mu and tau neutrinos. 
This should lead to a large muon number (number of mu neutrinos minus antineutrinos) or tau number of up to 10^^ 
for the hot protoneutron star 0|. This is accompanied by a nonzero chemical potential for mu and tau neutrinos. 

Accurate treatments of neutrino transport are now feasible based on the Boltzmann equation Therefore, it 

is important to systematically improve neutrino opacities. In section II, we catalog a number of possible opacity 
corrections. 

In section III, we present exact neutrino-nucleon cross sections for both charged and neutral currents. These are 
accurate to all orders in k/M where k is the neutrino energy and M the nucleon mass. We also present corrections 
from single nucleon form factors and strange quark contributions in the nucleon. In section IV we discuss recoil 
and weak magnetism corrections to mu and tau neutrino, energy and lepton number fluxes and spectra. Section V 
discusses corrections for electron antineutrinos. We also discuss possible changes in electron fraction Ye and heating 
rates. Finally, we conclude in section VI that these weak magnetism and recoil corrections should be incorporated in 
future supernova simulations. 



II. CORRECTIONS TO NEUTRINO OPACITIES 



Recently, accurate numerical algorithms have been developed to solve the Boltzmann equation for neutrino trans- 
port The high accuracy of these simulations now warrants systematic improvements in the neutrino opacities. 
Therefore, we discuss a number of opacity corrections. 

The corrections listed in table II can be divided into two main groups. Numbers 1 through 5 are classified as model 
independent because they can be calculated exactly and arc independent of the model used to describe dense matter. 
In contrast, corrections 6 through 13 are model dependent. Not only do these corrections depend on the model, but 
it is important that these corrections be consistent with the equation of state. 

Corrections 1, phase space, and 2, matrix element are perhaps the simplest. These are the only ones that are 
present in the limit of very low densities. They represent more accurate calculations of free space neutrino-nucleon 
scattering. By phase space we mean corrections from the Q value of the reaction such as the neutron-proton mass 
difference or because the outgoing neutrino momentum is slightly different from the incoming momentum. 

In this paper we focus on number 2, corrections to the matrix element. Recoil corrections, 2a, arise because the 
nucleon is not infinity heavy and recoils slightly from the neutrino. Weak magnetism, 2b, arise from the parity 
violating interference between the weak magnetic moment of a nucleon and its axial current. Weak magnetism is 
important because it has opposite sign for neutrinos and antineutrinos. Thus weak magnetism increases the opacity 
for neutrinos and decreases the opacity for antineutrinos. 

The first effects of the dense medium are 3, Pauli blocking of some outgoing nucleon or electron states, and 4, 
the Fermi or thermal motion of the initial nucleon or electron. These can be calculated exactly. The final model 
independent correction in table II is 5, coulomb corrections. If heavy nuclei are present, the opacity may be dominated 
by coherent elastic scattering from nuclei. However, these nuclei have strong coulomb interactions and form a strongly 
correlated classical liquid. The correlations between ions can greatly reduce the opacity for low energy neutrinos that 
have wavelengths comparable to the distance between ions. This correction is model independent because the coulomb 
interaction is known. For example, Ref. [pr] calculated the exact static structure factor of the ions with a simple 
Monte Carlo simulation. In principle, coulomb correlations are also present for neutrino electron scattering 

Corrections 6-13 depend on the model used for the strong interactions between nucleons. The mean field approxi- 
mation, 6, is perhaps the simplest way to treat the interactions. Nonrelativistic mean field models incorporate mean 
field effects in the density of states with an effective mass. Relativistic mean field models have both Lorentz scalar 
and vector mean fields and the scalar field changes the effective Dirac mass. 

It is important that the opacities be consistent with the model used for the equation of state (EOS). This can be 
achieved by using linear response theory [ p^ to calculate the response of the medium to a neutrino probe. The random 
phase approximation (RPA) or sum of ring diagrams provides the linear response of a mean field theory ground state. 
Therefore, if one uses a mean field EOS and calculates the opacities in RPA, where the RPA effective interaction is 
the same as that used for the EOS, then the opacities and EOS will be consistent. For example, the relativistic RPA 
opacities of Ref. are consistent with their mean field EOS. 
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We note that the nonrelativistic calculations of Burrows and Sawyer include corrections 1, 3-7 however they do 
not contain any matrix element correction, 2. Likewise the relativistic RPA calculations of |^ contain corrections 1, 
3-7 and 2a recoil but do not contain most of 2b weak magnetism. In the present paper we discuss weak magnetism 
corrections. In later work we will incorporate these into full nonrelativistic and relativistic RPA calculations. 

There are many NN correlations beyond a simple RPA approximation, 8, for example those responsible for NN ^ 
NNvv In addition there are meson exchange current corrections, 9. The Stony Brook group has done some 

work calculating corrections from hyperons or other strange hadron components, 10, p^ ]. More work could be done 
calculating opacities for meson condensates or other exotic phases of dense matter, 11. 

If superfluid or superconducting phases are present there will be pairing corrections to the opacity, 12. These can 



be important for the late time cooling of a neutron star 1 16| or perhaps if there is a color superconducting phase [n7^ 



Finally there could be important corrections to opacities in nonuniform matter, 13. For example, in the low density 
pasta phases present in a neutron star inner crust or for a nonuniform meson condensate [ p^ . 

All of the above assumed magnetic fields do not play important roles. If magnetic fields are important, they could 
substantially complicate the opacity, 14, see for example [^9). Magnetic fields might be important in jets or to break 
the spatial symmetry to explain neutron star kicks ■ 

To conclude this section, now that accurate neutrino transport is possible, it is important to systematically improve 
the neutrino opacities. One should include weak magnetism and other free space corrections along with the many 
density dependent efi^ects. Furthermore, these free space corrections can be included exactly and without any model 
dependence. 

III. NEUTRINO NUCLEON SCATTERING 

In this section we discuss recoil, weak magnetism, nucleon form factor, and strange quark corrections to the zeroth 
order cross section, Eq. (|l]). Then we present differential, transport and total cross sections that are exact to all 
orders in the neutrino energy over the nucleon mass k/M . 



A. Phase space 

Perhaps the simplest correction comes from the phase space of the scattered lepton because its energy k' is less 
than the incident energy k. If we use full relativistic kinematics for the nucleon and we assume it is originally at rest, 
then, 

fc' = fc/[l + e(l-a;)] EEryfc, (2) 
r^=[l + e{l-x)r\ (3) 

with X = COS0. 

We expect corrections to depend on the small parameter, 

e = k/M. (4) 

To the first order in e, 

k' K k[l - e{l ~ x)]. (5) 
If one uses the phase space with k' but still uses the matrix element for M ^ oo then Eq. (H) becomes. 



d<7ps G'^k''^ 



cl{l + x)+cl{3~x)], (6) 
-V', (7) 



dao 2 



dfl 

.|2|l-2e(l-.)l. (8) 

We note that Burrows and Sawyer's work [Q in the limit of low densities reduces to see Appendix I. 
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B. Recoil 



If one evaluates the matrix element to first order in e and includes the phase space corrections one has, 

We note that the full recoil correction in Eq. (^) is 50% larger than that for in Eq. (H). Equation (9) still 
neglects the parity violating interference between axial and vector currents. This interference is dominated by weak 
magnetism. 



C. Weak Magnetism 

The full weak vector current has Dirac (c„) and Pauli or tensor (F2) contributions, 



■ 2M 



(10) 



where q^^ — — k'^ is the momentum transfered to the nuclcon. The weak current is related by CVC (conserved 
vector current) to the electromagnetic current. Therefore the large anomalous magnetic moments of the proton and 
neutron give rise to a large F2. The couplings c„ and F2 are collected in table I. 
The cross section to the first order in e is, 



da 



dn 



± 



clil + x) + cl{3 - x) 



-3 e(l-a;) 



(11) 



with plus sign for v and minus sign for D. Weak magnetism increases v and decreases 9, while recoil decreases both 
v and P cross sections. Therefore weak magnetism and recoil corrections approximately cancel for v but add for 9. 
Note, Burrows et al. |^^, see also for example include weak magnetism for the charged currents but neglect it for 
neutral currents. Equation (11) is a good approximation for v. However it can become negative for large D energies. 
To cure this problem we include corrections to all orders in e. 



D. Exact cross section 

The exact cross section can be calculated, see also p2| , 

dh " l^'^ 1 tl + '^^ - ^(1 - ^)] +cl[l + v^+ v{l - x)] ± 2ca{c, + F2){1 ~ v') 

„2 



(1 - [F|(3 ~x)+ ic,F2{l - x)] I . (12) 



In general the F2 and CyF2 terms are small because they only enter at order e^. However, their inclusion is 
important to ensure a positive cross section at high energies. Note, Schinde considers the CaC„ term but neglected all 
F2 contributions p3|. 



E. Nucleon Form Factors 

The finite size of the nucleon can be included in Eq. (|l2|) by using appropriate form factors, 

F2 ^ F2{Q^), 

Ca^CaiQ^), (13) 
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that are functions of the four momentum transfer squared, 

= -ql = 2k^T^{l^x). (14) 

For example, for v^n —>■ e~P, 

Ca{Q^)^9a/{l + QVMir, (15) 

with ga ~ 1.26 and the axial mass Ma ~ 1 GeV. Form factors only enter at order e^. They reduce the cross 
section by about 10% at fc « lOOMeV, and even less at lower energies. However, for completeness we include simple 
parameterizations of the form factors in Appendix II. Because the corrections are so small, we will ignore the form 
factors in the rest of this paper. 



F. Total charged current cross section 



The opacity for Ve from charge current reactions depends on the total cross section for e p, 

da 



dfl 



dn. 



The zeroth order result is. 



2 1,2 



Co 



G^k 



while the correction factor from just phase space is. 



Rpsik) = cl 



i-^Ml + 2e) 



2e- 1 1 , „ , 

— H ^Zn(l + 2e) 

el + 2e 2e2 ^ ' 



l{cl + icl), 



R 



4 fcl + 5cl 



■ps 



3 \cl + 3c^ 



f e + 0(e2 



The exact result is, 



da 
dh 



dn = aoR{k), 
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R{k) ^{ci{l + Ae+ —e/j + 3c^, (^1 + -e j ±4 (c, + F^) CaC i^l + -e 
c.F,e' + (l + ^e) F|| / [{4 + 3cl) (1 + 2^)^ 

4 + 5cl+2 {Cy + F2) Ca" 



R{k) w 1 - 2 



For Ve using table I this is, 
while for h'e we have. 



3c2 



« 1 - 7.22e 



i?^ « 1 + l.Ole 



e + Oie/). 



(16) 

(17) 

(18) 
(19) 
(20) 

(21) 



(22) 
(23) 

(24) 
(25) 



We note the main effect of weak magnetism and recoil is to reduce the opacity by the large amount in Eq. (|24|). 
For a typical neutrino energy k ~< E >^ 20 MeV this is a 15% reduction. 

The full correction factor R{k), is plotted in Fig. ^ along with the lowest order forms in Eqs. ( p^ , p5| ). We see that 
the lowest order is fine for neutrinos but fails for antineutrinos above about 50 MeV. 
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G. Neutral current transport cross section 



For mu and tau neutrinos the scattering opacity depends on the transport cross section, 



(26) 



Note, strictly speaking this is true only in the limit k' = k nevertheless Eq. ( pq ) provides a good estimate when 
k' w k. 

The zeroth order transport cross section is, 



— 3(cS + 5c^) 



The transport cross section including only phase space corrections is. 



(27) 



(28) 



ii±^/n(l + 2e)-l 



2e+(2e2-e-l)ln(l + 2e) 



e3 (1 + 2e) 



(cl + 5cl) 



(29) 



Finally the exact transport cross section is 



2cg + Meg 



cr* = al R\k), 



e + 0(e2). 



(30) 



(31) 



R\k) 



e-1, , „ , 3 + 12e + 9e2-10e3 
ln{l + 2e) + 



2e3 
1 + e 



± (Ci, + F2) Ca 
1 



;ri(l + 2e) - 
1 



+F4 



lnil + 2e) + 



3e2(l + 2e)^ 

lOe^ + 27e2 + I8e + 3 

3e2(l + 2e)3 

, , „ , 2 + 10e+ f 

?n l + 2e ^ ^ .'3 

e(l + 2e)-^ 

2 /3 + 15e + 22e2 

3 V e(l + 2e)3 
' - 22e2 - 15e - 3 



ln{l + 2e) 



6e(l + 2e)3 
'2 



(cl + 5cl) 



(32) 



R* w 1 



-(3cg + 21c2)±8(c. +f2)Ca 
c2 + 5cl 



0(e2 



(33) 



The full correction, Eq. (|32[), is plotted in Fig. (g). Corrections for v ~ n are very close to those for v —p elastic 
scattering. We also plot the lowest order result, Eq. (|3^). This fails for antineutrinos at high energies. 
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H. Strange Quark Contributions 



Neutrino nuclcon elastic scattering is sensitive to possible strange quark contributions in the nuclcon. Strange 
quarks do not contribute to Cv{Q^ — 0) because the nucleon has no net strangeness. However both F2{Q^ = 0) and 
CaiQ^ = 0) can be modified. Parity violating electron scattering constrains F2 |2^. In this section we concentrate on 
Ca because F2 only contributes at order fc/M. 

Strange quark contributions are assumed to be isoscalar while the dominant contribution to Ca is isovector. There- 
fore, strange quarks are expected to increase the cross section for v — p scattering and decrease that for v ~ n. We 
write, 

Ca(Q' = 0) = i(±.g,-5:), (34) 

with ga = 1.26 and the plus sign for u — p and the minus sign for v — n scattering. The strange quark contribution 
5q is expected to be negative or zero. The best limit comes from a Brookhaven v ^p scattering experiment p^ . We 
assume, 

g:«-0.1±0.1. (35) 

A better measurement of from a laboratory experiment would be very useful. 

We make a simple estimate of strange quark effects starting from the zeroth order transport cross section in Eq. 
(p7|). We consider the neutral current opacity for a mixture of neutrons and protons of electron fraction Y^. For 
simplicity we set c„ = for v — p scattering. The ratio of the opacity with to that with 5^ = is, 

^^'^^ ^ ^9l + l-Y, ■ ^^^^ 

For = —0.2, S represents a 21% reduction in the opacity at Y^ = 0.1 or a 15% reduction at Y^, = 0.2. Thus strange 
quarks could reduce opacities by 10 to 20%. Note, strange quarks do not contributed to charged current interactions. 



IV. MUON AND TAU NEUTRINOS 

In this section, we use the weak magnetism and recoil corrections for transport cross sections, Eq. (|3^) to discuss 
muon and tau neutrino properties in supernovae. First, we examine muon and tau lepton number currents. Because 
neutrinos now have shorter mean free paths than antineutrinos, there will be net i^^ and densities and nonzero 
and Vr chemical potentials. We then calculate the increase in the and Vr energy flux because of weak magnetism 
and recoil. This should increase the neutrino luminosities. Finally we examine corrections to the emitted spectra of 
Vfi, Vr and v^^ Vr- We show that the P^j, Vt are expected to be hotter than i'^, Vr because of weak magnetism. 

For simplicity, we neglect charged muons. These could play some role at very high densities and or temperatures. 
Without muons, our results are identical for muon and tau neutrinos. Therefore, we write Vx in this section where x 
could be [i or t. 



A. Lepton number currents 

We start by considering high densities, well inside the neutrino sphere, where a simple diffusion approximation is 
valid. Earlier work calculated the lepton number current to first order in e = kjm. We extend this work to all 
orders. 

The transport mean free path for a Vx of energy fc is, 

with the zeroth order mean free path Aq evaluated at an arbitrary reference energy fcp, 

Ao= K(fco)p„]'\ (38) 
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and p„ is the neutron number density. Note, for simplicity we assume pure neutron matter, = 0. Because the 
correction factor i?^, Eq. (|3^), is so similar for v — p and v ~ n scattering our results are not expected to change much 
for Ye > 0. 

The lepton number current for in a simple diffusion approximation is, 

7(2^)3 3 ^l + exp[(fc-/i.J/r]- ^^^^ 
Here /i^^ is the neutrino chemical potential. For we have, 

^<'> ^ ^"FilW' <"»' 

This is now longer then A(fc) because of the different weak magnetism correction factors Rl < i?' . The lepton 
number current has the same form as Eq. with fi^,^ replaced by — /ij/^ , 

/_ =_ fjlLM^ I (41) 

J (2^)3 3 1 + exp [(fc + M.J/T] ■ ^ ' 

The total lepton number current is, 



+ cxp [(/c - /z,J/T] 
dfc 1 
i?Fl + exp[(fc + ^,J/r] 



(42) 



If J 7^ the lepton number of the star (number of Vx - number of v^) will rapidly change. This buildup in lepton 
number gives rise to a nonzero chemical potential /i^^ . 

Reference ^ argued that the system rapidly reaches steady state equilibrium where J = — Ju^ = 0. We 
numerically solve for the chemical potential so that, 

dk 1 dk 1 

(43) 



Ri{k) 1 + e('=-A'-x)/^ Jo Rl{k) 1 + eC^+A'-J/'T 

This is shown in Fig. |3| as a function of T . This chemical potential insures J — Q. 

If we expand Eq. ( [4^ ) to first order in [iv^jT and use Ri{k) and R\^{k) expanded to the first order in k/M one 
reproduces the lowest order chemical pot /io of ref. Q|, 

M..«Mo = -^^, (44) 

with 5c = 8(c,j + F^)ca/(c^ + 5c^) « 3.30 for v — n scattering. The lowest order /iq is seen in Fig. || to be an excellent 
approximation to the exact result even at T = 80 MeV. Note, this nonzero can lead to muon number (number of 
j/^ minus number of P^) or tau number for the protoneutron star as large as 10^'* ||l|. 



B. Energy Flux 



We calculate the energy flux carried by Vx by multiplying the integrand in Eq. (p3) by fc. 



d^k A(A:), 
(27r)3 ~Y' 



(45) 



We use Eq. (^) for th = aT"^ with a = (5c7r^/6M and add the Vx energy flux. We assume the temperature T 
has only a radial dependence so that. 



7 — 7B I tE _ ^o^o 



Gtt^ dr 



r dk f , e \ f 



^aT-k/T ^ 2 



Rh^T^ 



(46) 
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In the absence of weak magnetism and recoil corrections, iii,^ = and the energy flux becomes, 



•J R 



Gtt^ dr 



dk 



T2 



-k/T 



(47) 



In Fig. ^ we plot the ratio of the full energy current with weak magnetism and recoil to the zeroth order result 
JeI J%- We see that weak magnetism and recoil substantially enhances the energy flux. This should raise the neutrino 
luminosity of the protoneutron star. At a temperature of 30 MeV the enhancement is over 50%. 

Many supernova simulations combine and 9^; into one effective species. In this approximation /i^^ = 0. Therefore, 
we evaluate Je in Eq. ( ^ ) while setting a — 0. This would correspond to using a (geometric) average for the correction 
factor, 



<i?>-=- 



(48) 



for both i/x and Vx- We see that Jsia ~ ^)/Je slightly larger then Je/Je- This is because the chemical potential 
somewhat reduces the number of 9 and hence their energy current. Even so, J_B(a = 0) is a much better approximation 
then J'^. 



C. Spectrum 



The above lepton number and energy fluxes are based on a diffusion approximation. This is valid inside the neutrino 
sphere where A is much less than the size of the system. We now discuss weak magnetism and recoil corrections to the 
emitted Vx and Vx spectra. We will use a simple Monte Carlo model of Raffelt to estimate the emitted spectra. 
We emphasize that this simple model should be checked against full simulations that include weak magnetism and 
recoil. Nevertheless, we expect the model to provide a first orientation. 

Raffelt psi discusses the formation of spectra. At high densities, reactions such as e+e^ <-> Vx^x, VxC VxC^ or 
NN NNvxVx keep Vx in chemical and thermal equilibrium with the matter. We define the energy sphere as the 
approximate location where these reactions, which all have small cross sections, become too slow to maintain thermal 
equilibrium. 

Next the Vx propagate through a scattering atmosphere where neutrinos diffuse because of Nvx Nvx which has 
a much larger cross section. However the small energy transfer in nuclcon scattering is assumed to be too small to 
maintain thermal equilibrium. Finally, the neutrinos escape from the neutrino sphere. 

Note, the average spectrum for neutrinos and antineutrinos may be modified somewhat by nucleon-nucleon brem- 
strahlung in the scattering atmosphere pst . However, we expect weak magnetism in nucleon scattering to provide an 
estimate for the difference between the spectrum of neutrinos and antineutrinos. This should be explicitly checked in 
future work that includes weak magnetism for nucleon-nucleon bremstrahlung. 

We model the energy sphere as a black body with temperature Tes ~ 12 MeV. This high temperature is chosen so 
that the final temperature of the emitted Vx agrees with detailed simulations. Alternatively, Tes is the hot matter 
temperature well inside the neutrino sphere at densities high enough (-- lQ^^g{cm^) so that the Vx are thermalized. 
We assume the same Tes for both Vx and Vx because they are produced in pairs. 

We next assume the Vx propagate through a scattering atmosphere of optical depth. 



T{k) 



Toe 



Rl{k) 



(49) 



Here Tq is the thermally averaged optical depth in the absences of recoil and weak magnetism corrections. Raffelt 
assumes a Boltzmann distribution for the initial thermalized neutrinos for which < fc^ >= \2T^g. Note, the optical 
depth is the thickness of the scattering atmosphere in units of the energy dependent mean free path. A value of tq ^ 
30 corresponds to one typical simulation . 

Raffelt calculates the survival probability S{T{k)) for a Vx of energy k to make it through the atmosphere and 
escape. Otherwise it is assumed to be scattered back to the energy sphere and absorbed. A fit to ref. |^5| Monte 
Carlo results is. 



S{t) = 





-1 r 











1 - 



0.033 



1 + 1.5(4 + -17)2 + 0.5(4 + -32)6 



(50) 
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with (.r = logj^Q r. Therefore the final emitted spectrum is, 

UM = k\-''/^-^S{T{k)) (51) 

For antineutrinos the optical depth is, 

= TT^i?* (fc), (52) 

and the i>x spectrum is, 

« iVofc2e-'=/^--5(f(fc)) . (53) 

The normalization Nq is chosen so that dkf^^ = dkf^^ . This corresponds to steady state equilibrium and no 
net change in the lepton number. Note, Nq is related to the chemical potential /ii/^ . 

The i^x and the Dx spectra are shown in Fig. 5. Both of these curves differ from a Boltzmann distribution. 
Nevertheless the average energy can be characterized by a spectral temperature, 

T=<k>/3. (54) 

Table II collects T values. We find that T for i>x is about 7% larger than that for i/x almost independent of the choice 
of tq. Thus weak magnetism insures that the Dx spectrum is hotter than that for i/x- At high energies Rl{k) is small. 
Therefore fp{k) becomes significantly larger than f^^{k). This is shown in Fig. 5 as an increasing ratio of fp^ik) to 
fv^(k). Thus the high energy tail in the spectrum is expected to be antineutrino rich. 

The difference between the Vx and the Vx spectra in Fig. ^ is a macroscopic manifestation of charge conjugation, 
C, violation in the standard model. If the weak interactions had conserved C then i?^ — R^. In principle, this 
difference is directly observable. However, it may be very difficult to distinguish a detected Vx from a Dx- Instead, a 
superposition of the Vx and the Dx spectra may be more easily observable. Weak magnetism, by separating the Vx and 
the Vx spectra and by increasing the high energy tail of the Dx may lead to an observable broadening of the combined 
fx and Ux spectrum. 



V. ELECTRON ANTINEUTRINOS 



The reduction in charged current interactions from weak magnetism and recoil can change the emitted i>e 
spectrum, decrease the electron number current and increase the De energy flux. These changes, in turn, impact the 
Ye in the neutrino driven wind which is a possible site for r— process nucleosynthesis [|[^]. Finally these changes could 
be important for analyzing the detected Ve events from SN1987A or a future galactic supernova. 



A. Lepton number current 



The increased mean free path from recoil and weak magnetism will increase the De number current or decrease 
the total electron lepton number current Jg = Jj^^ — Jv^ ■ This can increase the density in the star. For simplicity, 
consider a situation where is near zero in a simulation without recoil or weak magnetism. Then in steady state 
equilibrium, we expect recoil and weak magnetism to lead to a /j,^^ of order the from section III Eq. (Q). This 
chemical potential will lead to an increase in the electron fraction because, in (3 equilibrium, /ig + = /ip + /ijy^. This 
will lead to a small change in the electron fraction of order. 



1 



(55) 



Here /x^ is the electron chemical potential and Y^ the electron fraction without weak magnetism and Yg is the new 
electron fraction. Using Eq. (|44|) for at T = 10 MeV, /x° = 15 MeV one has. 



Y,^Y^{1 + .12) 



(56) 



10 



an increase of 10%. Note, Eqs. (^,^^ slightly overestimate the change in because the charged current opacity for 
Pe increases with Yg. Nevertheless, weak magnetism could lead to a modest increase in Yg. 

In principle, weak magnetism will increase the lepton number diffusion time scale. However this time scale is very 
sensitive to temperature. We also expect weak magnetism to increase the Pe energy flux in a similar way to Fig. 4. 
With weak magnetism, the protoneutron star should cool faster and the lower temperatures will decrease the lepton 
number diffusion time. These contrary effects should be investigated in a full simulation. 



B. Spectrum of Vs 



We consider the following simple model [g5[ to estimate the change in the De spectrum with weak magnetism. Let 
the Dg neutrino sphere be at a temperature Tq without weak magnetism. For simplicity we assume a Boltzmann 
spectrum e~^/^°. 

Now we include weak magnetism and recoil described by Rc,{k). This reduction in the opacity will shift the neutrino 
sphere inwards to higher densities and temperatures. We define a new energy dependent temperature T{k) to include 
this shift, 



T{k) 



1 



Ru{k) 



Here u is the ratio of temperature to density gradients, 



dlnT 
dlnr 



dlnp 



rflnr 



To. 



1 



(57) 



(58) 



for the protoneutron star at the neutrino sphere. The extra -1 in the denominator follows because the optical depth 
involves a path integral from the neutrino sphere to infinity ||2^. Realistic values of u could be ~ 0.25 to 0.35 |p5| . 
The energy dependent temperature T{k) defines the emitted spectrum, 

/,^(fc)«/fc2g-fe/TW (59) 

If we expand Ri^{k) « 1 — Sck/M, with Sc = 4:{cy + F2)ca/ (c^ + 3c^) + 3 « 7.22 and expand /p^ to first order in fc/M, 

uSck"^ 



1 



MTo 



-k/To 



We define a spectral temperature Tp^ as, 

=< fc > /3 « To [1 + 8u6cTo/M] , 

For u = 0.3 and To = 5 MeV one has, 

Tp^ « 1.090 To. 



(60) 



(61) 



(62) 



Thus weak magnetism and recoil can increase the emitted De average energy by of order 10%. It is important to check 
our simple model with full simulations. 

Figure 6 shows the full spectrum, Eq. (|60|). We see that weak magnetism and recoil shift the strength to higher 
energies and increase the high energy tail significantly. This is shown by a large ratio of the spectrum with weak 
magnetism to k'^e~^/'^° at high energies. 



C. Conditions outside the gain radius 

We now discuss some possible implications of weak magnetism on neutrino heating, electron fraction and neutrino 
cooling in the low density region outside the gain radius. The gain radius is the point were cooling from neutrino 
emission balances heating from neutrino absorption and scattering. The net amount of heating outside the gain radius 
could be important for the success of the explosion. 
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The absorption of outside the gain radius will be reduced by the smaller cross section from weak magnetism. 
However the cross sections grows with the square of the energy so the new absorption rate may be proportional to, 



Rate - < -R^^elfc) > • (63) 

The first factor could increase the rate by about 20% while the second factor < Rp^ >~ 0.8 could decrease the rate 
by a similar amount. Therefore, the net rate of De absorption might be little changed. Thus the electron fraction, 
which depends on the absorption rate, should not change greatly. However, the DeS absorbed have a higher energy so 
the net heating should increase by about 10% because of the higher average energy in Eq. (|6^). 

This assumes the luminosity of the changed only because of the change in T^^ . If the luminosities of v,, and Df, 
increased further because of the increased energy fluxes in Fig. 4 then the total heating could be even larger. 

Finally we discuss cooling from positron capture outside the gain radius, 

e'^ + p ^ n + De- (64) 
The rate for Eq. (|6^) will be reduced by weak magnetism, 

-k/T 



poo 

Rate OC / dkk^kr' Ri,^{k)e~ 
Jo 



dkk^ [l-5r-\ e-'^/^, 



m ^ 

^5!r»(i-5^^). 

At T = 2 MeV the factor in parenthesis reduces the cooling rate by 9%. Thus weak magnetism can change the Ve 
luminosity to increase the heating while, at the same time, reducing the cross section to decrease the cooling. These 
results should be checked with full simulations. 



VI. CONCLUSION 



In this paper we examine recoil and weak magnetism corrections to z/-nucleon interactions. These are important 
because they are present at all densities, even at the relatively low densities near the neutrino sphere. Furthermore, 
the corrections are model independent. We calculate simple, exact correction factors Eqs. ( [T^j2^ , |3^ ) to include recoil 
and weak magnetism in supernova simulations. 

Perhaps the most important effect of weak magnetism and recoil is to increase energy fluxes of both and Vf. 
antineutrinos. We calculate, in a diffusion approximation, an increase of order 15% in the total energy flux for 
temperatures near 10 MeV. This should raise the neutrino luminosity. 

Weak magnetism and recoil will also change the emitted spectrum of and Ve- We estimate that will be 
emitted about 7% hotter than because have longer mean free paths. Likewise weak magnetism may increase 
the temperature by of order 10%. This increase in temperature coupled with the increase in neutrino luminosity 
should overcome the reduced absorption cross section and increase the heating in the low density region outside of 
the neutrino sphere. This, in turn, could be important for the success of an explosion. 

We find large corrections. However, supernova simulations are very complicated with many degrees of feedback. 
Therefore, it is important to check our results with a full simulation that includes Boltzmann neutrino transport and 
weak magnetism corrections. 
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Table I Coupling Constants 
Reaction c„ Ca F2 

vp vp 1 - 2 sin^ K. 0.035 S'„/2 « 0.630 - H„) - 2 sin'' ^^/Xp w 1.019 

i^n ^ — i — .9a/2 « -0.630 — ^(/Up — Mn) — 2sin^ ^-iuMn ~ --963 

ZlZ^n}^ ffa« 1.260 Mp-Mn« 3.706 

Here ga « 1.260, sin^ 6^ « 0.2325, /ip = 1.793 and /i„ = - 1.913. 



Table II Corrections to v Opacities 
Correction 



1. 




Phase space 


2. 




Matrix clement 




a. 


recoil 




b. 


weak magnetism 




c. 


form factors 




d. 


strange quarks 


3. 




Pauli blocking 


4. 




Fcrmi/thcrmal motion of initial nucleons 


5. 




Coulomb interactions 


6. 




Mean field effects 


7. 




NN Correlations in RPA 


8. 




NN Correlations beyond RPA 


9. 




Meson exchange currents 


10. 




Other components such as hyporons 


11. 




Other phases such as meson condensates or quark matter 


12. 




Corrections from superfluid/ superconductor pairing 


13. 




Nonuniform matter 


14. 




Magnetic field effects 



Table III Spectral temperature T =< k > /2> 
TO r., (MeV) Tp, (MeV) T^JT,^ 



3 8.75 


9.24 


1.056 


10 7.48 


7.96 


1.064 


30 6.49 


6.92 


1.066 


100 5.64 


6.04 


1.071 


300 5.09 


5.46 


1.073 
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APPENDIX A: LOW DENSITY LIMIT OF BURROWS AND SAWYER 



In this appendix we discuss the low density limit of Burrows and Sawyer's calculations and show that it reduces 
to our phase space result, daps/dft, in Eq. (^). Burrows and Sawyer start by considering a pure vector interaction. 
Their Eq. (2) for the differential rate of neutrino scattering k ^ k' is, 



dujdcosB 47r2 



(Al) 



Here the momentum transfered to the neutrons is q = k — k' and the energy transfered is lu = k — k' . The neutrino 
distribution is fu{k) and A''" is the neutrino tensor, Eq. (3) of ref. The dynamic structure function S'(g, w) for 
the nuetrons is. 



d^p 



f{p)[l - /(|p + q\)]27T6{uj + ep- ep+,), 



(A2) 



where f{p) is the neutron distribution function and ep the energy of a neutron of momentum p. 

One can recover our Eq. by, (a) neglecting the final state neutrino pauli blocking, (b) integrating Eq. (Al) over 
u), (c) assuming that the density is low enough so that 1 — /(p + q) « 1 in Eq. (A2), 



dujS{q, to) 



d^p 
(2^ 



/(p)27r. 



(A3) 



= 27rp„. 



(A4) 



Here p„ is the neutron density. Finally (d) we assume, in the limit of low density, that the lu depdendence of S{q, lu) is 
sharply peaked near, lu — k — k' , with k' given by our Eq. (^. Adding a simillar expression for the axial interactions 
gives, 



dT , ^ daps 



duudcosO 



dfl ' 



(A5) 



with dups/d^l given by our Eq. 



APPENDIX B: SINGLE NUCLEON FORM FACTORS 

In this appendix we collect simple parameterizations of the single nuclcon form factors for use in Eq. (|l|). We 
write the weak form factors in terms of the electromagnetic Dirac, Fl^ and Pauli, i^jj form factors for i = p or n. 

1) Reaction vp vp 



2 sin^ I 



iv) 



(«) 



(Bl) 



= |i (1 + 3.53r)-2 



F 



2= ^--2sm'e^\F^'" --F^ 



(B2) 
(B3) 



2) Reaction to 



c.^ I - -2sin^0^lF(")-^F(^) 



(B4) 
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Reaction Vf,n e p 



Here, see Eq. (14) 



with 



Ca = ffa (1 + 3.53t)" 



= F^P^ - 



F^P^ = [l + r(l + Ap)]G/(l + r) 



F^P^ ^\,G/{l + r) 



f(") « rA„ (1 - 77) G/(1 + t) 



and finally, 



F^"' = A„(1 + t,7)G/(1 + t), 



Ap = 1.793, A„ = -1.913, 



ry= (1 + 5.6t) \ 



G = (1 + 4.97t)" 



(B5) 
(B6) 

(B7) 
(B8) 
(B9) 

(BIO) 
(Bll) 
(B12) 
(B13) 
(B14) 

(B15) 
(B16) 

(B17) 



These parameterizations are from Ref. p6[ 
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FIG. 1. Factor R{k) that corrects the total charged current cross section for weak magnetism and recoil, Eq. (|22[), versus 
neutrino energy k. Solid line is for neutrinos, dashed line is for antineutrinos and the dot-dashed line includes only phase space 
corrections, Eq. (^^. Also shown as dotted lines are the corresponding lowest order correction factors, Eqs. (|2^) and (p3[). 



17 




k (MeV) 

FIG. 2. Factor R^{k), that corrects transport cross sections for weak magnetism and recoil, Eq. (^), versus neutrino energy 
k. The upper thick solid line is for vn — > vn while the upper thin solid line is for up i/p. The corresponding lower solid 
curves are for antineutrino scattering while the dotted curves give lowest order results, Eq. (|33| ) 



18 




T (MeV) 

FIG. 3. Muon or Tau neutrino chemical potential over temperature /i^^/T versus T for matter in steady state equilibrium. 
The solid line is the full result from the solution to Eq. (^3|) while the dashed line is correct to lowest order in k/M and /^/T, 
Eq. (@) 
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10 20 30 40 50 

T (MeV) 

FIG. 4. Ratio of the full energy flux of Vx and Ux neutrinos including weak magnetism and recoil Je to the energy flux 
without weak magnetism and recoil J% versus temperature. The solid line is the full result Eq. (Uq) over Eq. @. The dashed 
line neglects the neutrino chemical potential in Eq. and corresponds to = a = 0. 
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0.04 



2 




k (MeV) 



FIG. 5. Sp ectrum of muon or tau antineutrinos (solid curve) and neutrinos Vx (dashed curve) versus neutrino energy k, 
see Eqs. (O) and (M). The ratio of the Vx to Ux spectrum is shown by the dot dashed curve using the right hand scale. 
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FIG. 6. Spectrum of electron antineutrinos v^. versus neutrino energy k. The solid curve includes weak magnetism and recoil, 
Eq. (js^), while these are neglected in the dashed curve. The ratio of solid to dashed curves is plotted as the dot-dashed curve 
using the right hand scale. 
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